library(tidyverse)


crc230304 |>
  filter(POS == 161674008)

snpeff0304 <- read_delim('~/work/results/seekgene_crc/vcf/variant_annotation/crc230304/crc230304_snpEff.ann.vcf.gz', comment = '##')

snpeff0304 |>
  filter(ID == 'rs1050501') |>
  pull(crc230304)

snpeff0309 <- read_delim('~/work/results/seekgene_crc/vcf/variant_annotation/crc230309/crc230309_snpEff.ann.vcf.gz', comment = '##')

snpeff0309 |>
  filter(ID == 'rs1050501') |>
  pull(crc230309)

snpeff1125 <- read_delim('~/work/nextflowing/rnavar-crc221125/variant_annotation/crc221125/crc221125_snpEff.ann.vcf.gz', comment = '##')

snpeff1125 |>
  filter(ID == 'rs1050501') |>
  pull(last_col())
